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Abstract 

o 

OA ' In this work, we are interested in the controllabiHty of Vlasov-Poisson systems in the 

presence of an external force field (namely a bounded force field or a magnetic field), by 
' means of a local interior control. We are able to extend the results of [7], where the only 

present force was the self-consistent electric field. 

in 

1 Introduction and main results 

■ We consider the controllability of the Vlasov-Poisson system in the periodic domain T" (where n 

is the space dimension), which describes the evolution of a population of electrons in a neutral- 
izing background of fixed ions, under the influence of a self-generated electric field. The control 
questions are addressed by means of an interior control located in an open set co of the domain, 
which is a priori arbitrary. We assume in this paper that the charged particles evolve with 
the influence of an additional fixed external force, denoted by F{t, x, v) (at least with Lipschitz 
regularity and a sublinear growth at infinity in velocity). The equations read: 

00 ■ dtf + v.V.,f + F{t,x,v).V,,f + = 1^G, xGT", E (1.1 

p 
o 



A,$ = / fdv- fdvdx, (1.2) 

f\t=o = /o. (1.3) 



In these equations, f{t, x, v) is the so-called distribution function, which describes the density of 
particles at time t S M+, at position x G T" and velocity v £ M". The initial density distribution 
^ . fo(,x, v) is a non-negative integrable function. The right-hand side of the transport equation lujG 

I is a source term describing emission and absorption of particles, supported in uj. Moreover, to 

preserve global neutrality, G has to satisfy the following constraint: 

VteM+, / l^Gdvdx = 0. 

We normalize here the torus so that its Lebesgue measure is 1. 

The controllability problem is the following. Let fi{x,v) be another non- negative integrable 
function satisfying /i > and 

fidvdx = / fodvdx, 
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and let T > a fixed time. The question is: is it possible to find a control G such that: 

f{T,x,v) = h{x,v). (1.4) 

When the only acting force is the self-consistent electric field (that is when F = 0), the first 
author provided in [7] some positive answers to the question. More specifically, two kinds of 
results were obtained: first local controllability (which means that /o and /i are small in some 
weighted norm) were obtained in two dimensions, for an arbitrary control zone lo. Global 
controllability results (without restriction on the size of /q and /i) in any dimension was also 
obtained, provided that the control zone uj contains the image of a hyperplane of M" by the 
canonical surjection (which is called a hyperplane of the torus in [?])• The proofs of these results 
relied on the nice geometry of free transport in the torus: we shall recall their principle in a 
subsequent paragraph. 

When one considers a non-trivial external force the underlying dynamical system is more 
complicated; thus the characteristics can have a complex geometry, making the generalization 
not straightforward from the case F = 0. 

In this paper, we are able to extend results of for the two following classes of force fields: 

• The case of bounded force fields F G L'^Wx]^ . 

• In two dimensions, the case of Lorentz forces for magnetic fields with a fixed direction 
F(x,v) = b{x){v2, —vi) with b satisfying a certain geometric condition (which will be 
precisely described later). 

As we will see later on, the treatment of these two cases are rather different (in particular 
for what concerns high velocities) and involve different strategies. As a matter of fact, we were 
not able to find a general strategy which would allow to treat all forces F which are Lipschitz 
with a sublinear growth at infinity in velocity. 

Let us now briefly review the existing results on the Cauchy theory for the Vlasov-Poisson 
equation posed in the whole space R" or in the torus T". In this work, we will only focus on 
strong solutions (at least with a regularity in all variables); in the case where F = 0, the first 
results for such solutions are due Ukai and Okabe |10| who have proved global in time existence 
in two dimensions and local in time existence in three dimensions, in the whole space setting. 
One can readily check that the proof is the same for the torus case. In three dimensions, in 
the whole space, global in time results were proved independently by Pfaffelmoser [9] and Lions 
and Perthame |^ . The results of Pfaffelmoser were adapted to the torus case by Batt and Rein 
[3]. Concerning global weak solutions, the main result is due to Arsenev |T]. One can observe 
that all these results can be easily adapted to incorporate an additional external force F (with 
F satisfying the previous regularity assumptions). 

We will only rely on the construction due to Ukai and Okabe in the following. We are now 
in position to precisely state the main results proved in this paper. 

1.1 Results in the bounded external field case 

We first consider the case where F G L^Wx;v. In this case, we are able to exactly extend those 
for F = 0, that are a local and a global controllability results. The local result concerns only 
the dimension n = 2, but is valid for any control zone lo. On the contrary, the global result is 
valid for any n, but requires a stronger geometric assumption on the control zone u. 
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Theorem 1.1 (Local result). Let n = 2. Let F{t,x,v) e LfW^;^. Let > 2 and T > 0. 
There exist k,k'>0 small enough such that the following holds. Let /o and fi be two functions 
in C1(T2 X M2) n T^1'°^(T2 X R"^), satisfying the condition that for any {x,v) G x and 
iG {0,1}, 

r \f,{x,v)\<Kii + \v\r^-\ 

\ \yxfi\ + \VJi\<K'{l + \v\)-\ 

and 

[ fo= [ fi. (1.6) 

Then there exis ts a control G £ C°([0,T] x x R^), such that the solution of n~l\) - [r^) and 
I exists, is unique, and satisfies ([T]7P- 

Theorem 1.2 (Global result). Let 7 > n and k,k' > 0. Suppose that the regular open set uj 
contains the image of a hyperplane in M" hy the canonical surjection, supposed to be closed. Let 
/o and fi be two functions in C^(T" x M"), satisfying the conditions 

r mx,v)\<K{i + \v\r^-^ 

\ \^xfi\ + \yvfi\<n'{l + \v\)-\ 

and U.6\) . Then there exis ts a control G G C°([0,r] x T" x W), such that the solution of 
lll.l\) - [T^) and exists, is unique, and satisfies l{1.4\ )- 

1.2 Results in the magnetic field case 

Let us now state our result when F represents an external magnetic field. For all results dealing 
with this case, we will systematically assume that the space dimension n = 2. First, let us 
explain the physical meaning of the system under consideration. In the physical space M'^, let 
(61,62,63) a fixed orthonormal base. We consider the stationary magnetic field B, with fixed 
direction 63: 

B{x) = b{x)e3, 

where 6 is a Lipschitz function on T'^. Since B has to satisfy the divergence free condition, this 
implies that b only depends on xi and X2- The associated Lorentz force writes: 

F = V AB{x) = b{x)v^, 

denoting = {v2, — fi,0). We then restrict to distribution functions which do not depend on 
X3 and f3, so that we can restrict the study of the dynamics to the bidimensional plane (ei, 62). 
For the sake of readability, we rewrite the Vlasov-Poisson system that we study: 

dtf + v.VJ + b{x)v^.V,f + V,^.V,f = l^G, xGT^ vGR" (1.8) 
A^$= / fdv- [ fdvdx, (1.9) 



f\t=o = fo- (1.10) 
We now precisely state the geometric assumption we have to make on b. 

• Fixed sign. We assume that b has a fixed (say non-negative) sign. 
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Geometric control condition. We assume that there exists K a compact set of on 
which 6 > and which satisfies the geometric control condition: 



For any x G and any direction e G S 
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there exists y G M"*" such that x + ye ^ K. (1.11) 



One can notice that the geometric control condition corresponds to the geometric control 
condition of Bardos, Lebeau and Rauch [2] for the controllabihty of the wave equation. Let us 
underhne however that here this condition concerns the magnetic field only, and not the control 
zone u. As we will see, this condition assures that the particles are sufficiently influenced by the 
magnetic field. 

Examples. Let us give some examples, where this geometric assumption is satisfied. 

1. The most simple example that one can have in mind is the case where b is positive on T^. 
Then taking K = 1"^, the geometric assumption is satisfied. Obviously, this includes the 
case where 6 is a positive constant. 

2. Assume that h is non-negative and has finite number of zeros xi,...,xn G T^. Then 
there is r small enoug h such that K = T\ uf^^ B{xi , r) is appropriate. One could also 
extend this consideration to the case where the zeros of b are given by a sequence (xi)jgN 
with a finite number of cluster points. 

3. We can consider some b which is identically equal to in a large set of the torus, provided 
the existence of some K satisfying the geometric control condition. For instance, if we 
identify with [0, 1]^ with periodic conditions, a subset K containing ({0} x [0, 1]) U 
([0, 1] X {0}) satisfies the geometric assumption. 

With these particular magnetic fields, we are able to prove a local controllability result, 
which is similar to Theorem 11.11 (but we emphasize once again that the proofs will be rather 
different). 

Theorem 1.3. Let b satisfying the geometric assumption (jl.lip . Let 7 > 2 and T > 0. There 
exist k,k' > such that the following holds. Let fo and fi he two functions in C^iT"^ x M^) n 
1^1,00^^2 ^ ]^2-j^ satisfying the condition that for any {x, f;) G x and i G {0, 1}, 



Then there exists a control G G C''([0,T] x x M^) , such that the solution of \L8^ )- [T7§^) and 
il.lO\) exists, is unique, and satisfies f{T,x,v) = f\. 

1.3 Organization of the paper 

The paper is organized as follows: first, in Section [21 we will recall some considerations on the 
Vlasov-Poisson equation and will explain the general strategy of the proofs. Then, we prove 
Theorem 11.11 in Section [3] and Theorem 11.21 in Section |4l for what concerns the bounded external 
field case. Finally, in Section [Sj we prove Theorem 1 1.3 1 on the local controllability in the external 
magnetic field case. 




(1.12) 



and 




(1.13) 
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2 Strategy of the proofs 



2.1 Notations 

For T > 0, we denote Qt := [0,T] x T" x R", and Qt ■= [0,r] x T". For a domain we write 
also Clin), for I G N, for the set C\n) n W^^°°{n). All the same, C^+'^(fi) for a E (0, 1) stands 
for the set of functions with bounded cr-Holder /-th derivatives. Also, (7^''+'^ (f^r) (resp. 
C^'^^^ (Qt)), for / G N, o", fj' € [0, 1) is the set of continuous functions in (resp. Qt), which 
are with respect to x (resp. to {x,v)), and which /-th derivatives are all with respect to 
t and with respect to x (resp. to {x,v)). 

For X in T" and r > 0, we denote by B(x,r) the open ball with center x and radius r, and 
by S{x, r) the corresponding sphere. The radii will always be chosen small enough in order that 
S{x,r) does not intersect itself (that is r < 1/2 in the standard torus). 



2.2 The case F = 0, obstructions to controllability 

In this paragraph, we focus on the case F = 0, following [7]- Let us consider the linearized 
equation around the trivial state (/, $) = (0, 0). The linearized equation happens to be the free 
transport equation, which simply reads: 

dtf + v.VJ = Kg. 

By Duhamel's formula, we obtain the explicit representation for /: 



f{t,x,v) = fo{x-tv,v)+ / {lujg){s,x - (t - s)v,v)ds, 

Jo 

from which one can observe that there are two types of obstruction to controllability: 



(2.1) 



(small velocities) The second obstruction concerns the small velocities. The velocity of a 
particle can have a good direction, but if it is not high enough, then it will not be able to 
reach zone in the desired time, see Figure [TJ 

(large velocities, wrong direction) The first obstruction is of geometric control type as in 
[2] for what concerns the wave equation: if a particle has initially a wrong direction, then 
it will never reach the control zone, and thus we cannot influence its trajectory, see again 
Figure [TJ 



Figure 1: Obstructions for small and large velocities 

It follows that in general, the linearized equation fails to be controllable. 
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2.3 The return method 



In order to overcome these obstructions, the idea is to use the return method, which was intro- 
duced by Coron in [3] for the study of the stabihzation of finite-dimensional systems, and then 
used in the context of the control of PDEs by Coron in [Sj for the control of the two-dimensional 
Euler equation for perfect incompressible fluids. It has been used since in many different contexts 
of PDE control: we refer to the monograph of Coron [6] for several illustrations and references 
for this method. The principle is to build a reference solution (/, starting from (0, 0) and 
reaching (0, 0) in some fixed time, and around which the linearized equation enjoys nice control- 
lability properties. Such a construction can be delicate, and crucially depends on the structure 
of the studied equation. 

Here, the problem is more or less equivalent to find solutions (/, $) (starting from (0, 0) and 
reaching (0,0)) and such that the characteristics associated to V<& satisfy: 

VxeT", VvEM", 3te[0,T], X{t,0,x,v) Geo. (2.2) 

(As a matter of fact, the characteristics will not be quite associated to V$ inside the control 
zone.) 

When no exterior force is present, the existence of such a reference solution / was proved by 
the first author in in two dimensions, for an arbitrary control set lo. This is achieved using 
complex analysis tools by building harmonic potentials outside u, which allow to sufficiently 
influence the trajectories, so that the two previous obstructions are circumvented. This strategy 
distinguishes between high and low velocities, for which the relevant potentials are different. 

2.4 On the scahng properties of Vlasov-Poisson equations 

We notice that (|l.ip - ()1.2p is "invariant" by some change of scales. More precisely, when / is a 
solution of (HIID-dO]) in [0,r] x T'' x M", then for A / 0, the function 

f\t,x,v) ■.= \\\^-^f{Xt,x,v/\), (2.3) 

is stih a solution of ([ri])-([r2]), in [0,r/A] x T" x M" for the following potential 

ip^{t,x) := X^ip{Xt,x). (2.4) 

and the external force 

F^{t, X, v) := X^F{Xt, x,v/X). (2.5) 
The choice of some particular parameters A will be of great help for the controllability problem. 

The choice A = — 1. Using ()2.3p with A = — 1, we observe that in order to prove Theorems 11.11 
and 11.31 it is sufficient to prove the result for the case where /i = in [T"'\w] x M". Indeed, 
we observe that after imposing ()2.3p with A = — 1, the corresponding external field remains in 
the same class, that is, if F is bounded, then F^^~^ is still bounded (resp. if F corresponds to 
a magnetic field satisfying the geometric condition, then F^^~^ still corresponds to a magnetic 
field satisfying the fixed sign and the geometric conditions). 

Then one can follow the procedure that we detail below: 

- Take /q as initial value and (in (T"\tj) x M") as the final one, 

— Take {x,v) i— )• fi{x, —v) as initial value and again as the final one within the force field 
F{T -t,x,-v). 
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each in time T/3. We obtain two functions /o and /i. Now we may consider the function / 
partially defined in Qt by 

f{t, X, v) = fo{t, X, v), in [0, T/3] x T" x M", 
f{t,x,v) = 0, in [r/3,2T/3] x [T"V] x W\ 
f{t, X, v) = /i(r - t, X, -v) in [2r/3, T] x T" x W. 

Then we can complete in a regular manner / inside [r/3, 2T/3] x a; x M", taking care to preserve 
for any t the value of /u-nxK" /(^j v)dxdv. Finally we get a relevant solution /. For this reason, 
we will systematically assume that /i = in [T"\(x!] for all controllability results discussed in 
this work. 

The choice < A ^ 1: The choice of the parameters in such a range is useful to prove global 
controllability results. As in [7], it will help us in particular to prove Theorem 1 1.2 1 in the bounded 
external field case. The principle is that when A is chosen small enough then V^>^ has a small 
L°° small (and this is also the case for F"^), so that we can expect characteristics for to be 
close to those of some well chosen relevant reference solution. This will allow us to get rid of the 
smallness assumption on /q. Nevertheless in order to avoid concentration effects, we will need 
some assumptions on the characteristics associated to the reference solution. 

In the magnetic field case, we observe that F{x,v) = b{x)v'^ and thus F^{x,v) = Xb{x)v'^. 
For this reason, due to our treatment of high velocities for this case, this will not allow us to 
prove a global result. 

2.5 General strategy for external force fields F 

Following [7] the main steps for proving local controllability results will be: 

Step 1. Build a reference solution (/, of (jl.ip - p.2p with a certain control G, starting from 
(0,0) and arriving at (0,0), such that the characteristics associated to F — V<& satisfy ()2.2p . 

Step 2. Build a solution (/, $) close to (/, $), taking into account the initial condition (/o, <&o) 
and still arriving at (0, 0) (outside uj). This is achieved using a fixed point operator involving an 
absorption process in the control zone. This is where we use the smallness assumption on /q. 

The treatment of Step 2. will be quite similar to that in [7], although a bit more technical since 
we will have to take into account the geometry due to F. The main difference is the treatment 
of Step 1., for which we have to propose new ideas. The strategy is the following: 

Bounded force field. Our strategy relies on the fact that for short times, the dynamics with 
the external force F is well approximated by the dynamics with F = 0. We recall that in [7], 
the reference solution can be constructed for any time (which can be arbitrarily small) and any 
control zone in the torus. Thus, we use the construction in the case F = 0, for very short times 
and a small subset of the control zone uj, and using the approximation of the dynamics, this will 
give us a relevant reference solution. 

Magnetic field. The strategy in this case can be understood in the most simple case, that is 
when 6 is a positive constant. In this case, the characteristics associated to the magnetic field 
can be explicitly computed: these are circles, whose radius is proportional to the norm of the 
velocity (which is a conserved quantity). We make two crucial observations: 

- When the velocity is very large, the curvature of the circles are close to zero, and at least 
locally (that is for small times), the trajectory is well approximated by the straight lines 
of the free transport case. 
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— The magnetic field has "mixing features", in other words it makes the velocities of particles 
take every value of S^, which removes the above obstruction concerning high velocities. 
Hence, due to this effect, at high velocity, we do not need to create any additional force 
field to make the particles cross the control zone. 

This means that at high velocity any subset lo of the torus automatically satisfies the geometric 
condition (jl.lip for the caracteristics associated to the magnetic field. 

In the general case, the geometric condition on b allows us to make sure that the particles are 
sufficiently influenced by the magnetic field, so that the previous considerations will still hold. 

2.6 On the uniqueness of the solution 

In this paragraph, we briefly discuss the uniqueness question included in the above results. 

The flrst point is that, if we drop the uniqueness from the conclusions of the above theorems, 
we can replace the assumption 

WM + Wvfi\<^'{i + \v\)-\ 

by the weaker one 

\^xfi\ + |V„/i| < k'. 

This is easily seen when reading the proofs below. 

Hence the assumptions is of V/j belonging to some weighted space is only useful for the 
uniqueness issue. Let us explain how one can show uniqueness under this assumption. The 
main point is that in this case the solution described above satisfles 

\y .,vf{t,x,v)\ < C{h, h){l + \v\)-\ 

for all t. This follows from the construction described below, and from the estimates on V/ in 
the proof. Once these estimates are obtained, the proof of uniqueness is exactly the one of Ukai- 
Okabe. It consists in making the difference of two potential solutions; this difference satisfies a 
certain transport equation with source. Then one performs an n L°° estimate on the solution 
of this equation and uses a Gronwall argument. In our case, the source term disappears when 
we make this difference, so one can follow |10j without change. 
This gives the uniqueness among the solutions satisfying 

/GCi([0,r] xT^xM"), \f\ + \V^^,f{t,x,v)\<C{l + \v\)~'< and V(/. G L°°(0, T; 1^^'°° (¥")). 

3 Bounded external field case 

In this section, we prove Theorem 11.11 As already explained, the main difficulty is to build the 
reference solution. Then one can use the same absorption process, that was proposed in [7], and 
find a solution to the non- linear system by a similar fixed-point argument. 

3.1 Design of the reference solution for the bounded field case 

We begin with the construction of the reference solution. Accordingly to the previous strategy, 
we distinguish between high and low velocities. 

For the large velocities, we prove the following proposition: 
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Proposition 3.1. Let t > and H £ L°°((0,r) x T^;]^^). Given xq in and tq a small 
positive number, there exist tp G C°°([0,r] x T^;IR) and m > such that 

Aip = m [0,t] X [T2\;B(xo,ro/10)] (3.1) 

Supp(/9 C (0,r) X (3.2) 

and such that, if one consider the characteristics {X, V) associated to the force field H + V(/? 
then for all m > m: 

Vx G T^Vv G such that \v\ > m, 3t £ {t/3,2t/3), 



such that X{t,0,x,v) G B{xQ,rQ/4) and \V{t,0, x,v)\ > — . (3.3) 

Proof of Proposition I3.il In the case H = 0, this proposition was ah^eady proved in [TJ Propo- 
sition 1, p. 340]. We fix x'q = xq, Tq = ro/2. Applying this result for r = 1, we thus obtain 
the existence of (pi G C°°{[0, 1] x T'^;M) and m' G with compact support in time in (0, 1), 
satisfying: 

= in [0, 1] X [T\B{xo, ro/20)] , (3.4) 
Supp(^ C (0, 1) X T^, (3.5) 
and such that, if one consider the characteristics {X^ , V^) associated to the force field Vfi then: 

Vx G T'^,\fv G M"^, such that \v\ > m,3t e (1/4, 3/4), X^i, 0, x, v) G B{xo,ro/8). (3.6) 

Let r' < r to be fixed later. For this given r', we can construct (Pt-i by rescaling ipi as follows: 

(Pr'{t,x) := j^ipi(^^,x^ , (3.7) 

which corresponds to follow the characteristics with time ^. 

Now let us consider the shifted in time potential (p defined by: 



t-t' 



(p{t, x) = (pr> [t — , X ) . (3.8) 

We extend 99 by in (0, r) \ (^^, ^ 

We define the characteristics {X ,V) associated to the force field Vy?, which satisfy by con- 
struction: 

Vx G T'^',\/v G M'^, such that > m,3t e ( ^~^ , ^ t ),X{t,0,x,v) G B{xo,ro/8). (3.9) 



Let us now compare {X,V) and (X,V), which is associated to the force field H + on 
(0, r). By Taylor's formula we have: 

T — t' ~ t — t' /"* r ~ T — t' T — t' 

\X{t, — ^ — ,x,v) - X{t, — ^ — ,x,v)\ < 1^^ {t - s)^\\/ip{s,X{s, — - — ,x,v)) - Vip{s,X{s, — - — ,x,t;))| 



T — t' - T — t' 

+\H{s,X{s, — - — ,x,v),V{s, — - — ,x,v))\ 

(3.10)' 



ds. 
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By Gronwall lemma we deduce for t G i ■ 

\V{t, ^^,x, v) - V{t, v)\ < T'||i/||L^^_^e^"^'^"^°°((o>-)xT'') , 

The crucial point is now to observe that if described above satisfies: 



(3.11) 



|VVllL-((o,r)xT'i) = C" (^;^) asr'^O, 



as it can be seen from p.7p . 

Thus for t' small enough we infer that X{t, 0, x, v) meets B{xQ,rQ/4:) for some t G (^ ^^^ , ^~2^) 

(r/3, 2r/3), for all x and f, provided that \V{^^-^^,0,x,v)\ is large enough. This is ensured if 
\v\ > m is chosen large enough, thanks to the inequality: 

\V{^,0,x,v)\ > \v\ - ^\\HU^^^^. 

□ 

Remark 3.1. In this proof, this is crucial that H G L°°((0, r) x T^;R^). Thus this approach 
will fail for the magnetic field case. 

The above proposition shows that with a suitable electric potential, all particles having a 
sufficiently high velocity will eventually reach uj. The following proposition explains how one 
can accelerate all particles in order to make all the remaining ones also reach u. This will also 
rely on the construction in the case F = 0. 

Proposition 3.2. Let t > 0, M > and H e L°°((0,r) x T2;R2)_ g'j^g^ j2 ^ 

small positive number, there exists M > 0, £ € C°°{[0,t] x T^;R'^) and ip G C°°([0,r] x T^;R) 
satisfying 

£ = Vip in [0,r] x {T\B{xo,ro)), (3.12) 
Supp{£) C {0,t) xT"^, (3.13) 
Aip = m [0, r] x {T'^\B{xo,ro)), (3.14) 

such that if {X, V) are the characteristics corresponding the force 

I:=£ + H, (3.15) 

then 

y{x, w) G X B{0, M), V{t, 0, x, v)) G B{0, M) \ B{0, M + 1). (3.16) 
Proof of Proposition COi By [71 Lemma 3, p. 356], there exists 9 G C°°(T^;]R) such that 

A0 = O in T'^\B{xo,ro), 
\Ve{x)\ >0 in T^\B{xo,ro). 

From the second condition, one sees that Ind^^^.^, ,,p)(V0) = 0, so that '^0\j2\B{xo,ro) 
extended to as a smooth non-vanishing vector field, let us say W. Call A G C^((0, 1);M) a 
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nonnegative function with A = 1. We claim that for sufficiently small t' < t, and sufficiently 
large C > 0, 

Sit,x) ■■=^^{^) Wix), 

is convenient. Then all properties above but (|3.16p are clear. 

Call {X, V) the characteristics associated to £ only. We see that for all (x, v) £ T'^ x B{0, M) 
and t G [0,r'], 

|F(t, 0, X, v) - v\ < C||£:||oo, \X{t, 0, X, v)-x\< t'(C||£:||oo + M), 



so 



\V{t',0,x,v)-v + C£{x)\ < T'||f||,[r'(C||£:||oo + M)]. 



Noting that, due to the time support of £, V{t, 0, x, v) = V{t', 0, x, v) and using that \£\ > c > 
on T^, one sees that one can choose C and then r' such that 

V(3;, G X B{0, M), F(r, 0, x, v) eR^\ B{0, M + 2 + r||i?||oo). 

We now consider the characteristics {X, V) associated to £ + H and evaluate: 

\X{t,0,x,v) - X{t,0,x,v)\ < [ \V{s,0,x,v) -V{s,0,x,v)\ds 







\V{t,0,x,v) -V{t,0,x,v)\ < J (^\£{s,X{s,0,x,v)) - £{s,X{s,0,x,v))\ 

+ \H{t, X{s, 0, X, v),V{s, 0, X, v))\) ds 



<l|V<f ||L«>((0,r')xT2) / it- s)\V{s,^,X,v) -V{s,^,X,v)\ds + t\\H\\Lf 

Jo 

(3.17) 

By Gronwall's inequality: 

\V{t, 0, X, v) - V{t, 0, x,v)\< t\\H\\Lf^^yi\\^^\\ (3.18) 
We observe that we have: 

— ||V£:||ioo((o,,oxT2) = O(t') asr'^0. (3.19) 
Taking r' small enough, using t = r' in (|3.18p . and observing that 

\V{t,0,x,v) -V{t',0,x,v)\ < |r-r'|||i7||oo, 

allow us to prove our claim. The existence of M is a matter of compactness of x B{0,M + 
2 + t\\H\\^). □ 

Remark 3.2. We can observe that there is some "margin" in the previous proof, in the sense 
that if we only had 

^l|V<?llL-{(o,r')xT2) = C'(l) as t' 0, 

the proof would still follow. However, that (|3.19p holds will actually be crucial in the proof of 
the equivalent lemma in the magnetic field case, and this time this will be sharp. 
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The reference solution. Now we are able to define the reference solution. Consider xq in lo 
and ro a small positive number such that 

B{xo,2ro) C uj. 

We first define a reference potential Tp : [0, T] x — )■ M as follows. We apply Proposition 13.11 
with r = r/3, H = -F|[o,t/3]i we obtain Tpi and some > such that (|3.3p is satisfied. 
Let 

a = max (^^^^^, Croi^ + \\F\\oo + ll^illoo + ll^slloo)^ , (3.20) 
T 

Ml = max(mi,2a) + - (||V^i||oo + ||i^||oo) , Ma = max(m3,2a), M = max(Mi,M2). 

(3.21) 

Above Cro is a positive geometric constant depending only on rg, and which will be described 
later. 

We also use Proposition 13.11 again with r = T/3, H{t,x) = F{t + ^^x) for t G [0, T/3], we 
obtain Tp^ and some > such that ()3.3p is satisfied. Then we apply Proposition 13.21 with 
r = T/3, H{t,x) = F{t + j,2;) for t G [0,T/3], and M described above. We obtain £2 ^2 ^^'^^ 
some M. 

Finally we set: 

(^i(t,-) forte [0,f], 
^(i,-) = <; (?2(i-i,-) fort G [f,f 



(^3(t-f ,•) forte [f,T], 



and 



V(^i(t,-) fort G [0,^], 
£{t,-)={ £:2(t-f,.) forte [|,f], 

V(?3(t-f ,•) forte [f,T]. 

Let us now introduce /. Consider a function Z € C^(M";M) satisfying the following constraints 

' Z > in M", 
Supp 2: C -BiRn(0, 1), [2,22) 
I Z = l. 

We introduce / = /(t, x, v) as 

7(t, X, u) := 2:(u)A^(t, x). (3.23) 
Of course, / satisfies in [0,T] x x M^^ with source term 

G(t, X, v) := dj + v.Vj +{F + Vlp).Vj, (3.24) 

which is supported in [0,T] x B{xo,rQ) x M?. Up to an additive function of t, the function ip 
satisfies the equation (jl.2p corresponding to / (with /(O, •, •) =0). We denote 

p{t,x) := / f{t,x,v)dv = Aip{t,x). 
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3.2 Fixed point operator 

To prove Theorem II. H we construct directly the solution / starting at /o and reaching in 
\ cj at time T, provided that /q is suitably small. This is done by a fixed-point procedure. 
In this subsection, we describe the operator; in the next ones, we will find a solution to our 
controllability problem as a fixed point of this operator. 

Let e S (0, 1). We first define the domain of by 



Se-.-- 



< e, 



a. \\ luAa - f)dv\\csi^n^ 

b- 11(1 + \v\V{g - Dh^iQr) < Cl ll/o|lci(Tr2xM2) + 11(1 + ^1 )^/o IIcO(T2 xR2) 

C- Wg - f\\c^^2^Q^) < C2 ll/o|lci(T2xR2) + 11(1 + bl)'^/o|lcO(T2xR2) 

d. VtG[0,r], J.^2^^2g{t,x,v)dxdv = fj2y,^2foix,v)dxdv^, 



(3.25) 



with ci, C2 depending only on 7, T, u (and hence on {f,ip)) and but not on e. The indices 
5i < 82 in (0, 1) are fixed as follows 



-51 := 



7 — n 



and do 



7 



2(7 + 1) ' " 7 + 1 

For fixed ci and C2 large enough depending only on {f,{p), and /q small enough, one has 



(3.26) 



fodvdx 



and consequently, in this case /q + / € 5e, so 5^ 7^ 0. From now, this is systematically supposed 
to be the case. 

Now we introduce the following subsets of 5(xo,ro) x E?: 



1 



1 



7 := < {x,v) G 5(xo,ro) x M / k| > 7: and v.i'{x) < — — \v 



10' 



{x,v) G S{xo,ro) X / |u| > 1 and v.u{x) < —-\v\ 



(x,v) G 5(xo,ro) X / Iwl > 2 and v.iyix) < — If 

5 

7+ := {(x,t;) G ^(xcro) x / ^;.;^(2;) > o} , 



(3.27) 

(3.28) 

(3.29) 
(3.30) 



where i^ix) stands for the unit outward normal to the sphere S{xQ,rQ) at point x. It can be 
easily seen that 

dist([5(xo,ro) xM2]\^2-.^3-) 



We introduce a C°° fl regular function U : -^(xo, ro) x — ;> I 

< [/ < 1, 

C/ = l in [5(xo,ro) xM2]\^2-^ 
J7 = in 7^-. 



i, satisfying 



(3.31) 
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We also introduce a function T : 



M+, of class C°°, such that 





T' 




r47r 1 


in 




U 






"■48. 


48 ' 



and T = 1 in 



T 23r 



24' 24 



(3.32) 



Now, given g ^ S^, we associate ^p^ on [0,T] x by 



A(p3{t,x) = Jj^„g{t,x,v)dv - Jj„^^„g{t,x,v)dvdx in [0,r] x T", 
Jj„ ip^{t, x) dx = in [0, T]. 



Then, we define V{g) := / to be the solution of the following system 
' f{0,x,v) = foonJ^ xR\ 

< dtf + v.V^f + {F + VipS + g -VTp).V,f = m [0, T] x [(T« x M«)\7-], 
_ f{t,x,v) = [l-T{t)]f{t~,x,v) + Tit)U{x,v)f{t~,x,v) on [0,T] X7-. 



(3.33) 



(3.34) 



To explain the last equation, we introduce the characteristics (X, V) associated to the force 
field F + V(/?^ + £ — VTp. In the previous writing, f(t~,x,v) is the limit value of / on the 
characteristic {X,V){s,t,x,v) as the time s goes to t~ . (For times before t, but close to t, the 
corresponding characteristic is not in 7~.) When the characteristics {X,V) meet 7" at time t, 
then the value of / at time t"*" is fixed according to the last equation in ()3.34p . One can see the 
function T{t)U{x,v) as an opacity factor which varies according to time and to the incidence 
of the characteristic on S{xQ,rQ). In this process a part of / is absorbed on 7", which varies 
from the totality of / to no absorption according to the angle of incidence, the modulus of the 
velocity and the time. 

The set of times when a characteristic meets -y~ is discrete. Indeed, if {X, V){t,0,x,v) € 7" 
and {X,V){t' ,0,x,v) S 7", then there exists s E it,t') for which {X,V){s,0,x,v) E 7"^. The 
conclusion follows from dist(7"'", 7") > 0. 

We now consider a continuous linear extension operator vf : C^{T'^\B{xo, 2ro); M) — )• C^(T'^; M) 
and which has the property that each C"-regular function is continuously mapped to a C"- 
regular function, for any a E [0, 1]. 

From this operator, we deduce a new one vr : C°((T"\5(xo, 2ro)) x M") CO(T" x M") 
according to the rule: 

{7rf){x,v):=[Wf{;v)]{x). (3.35) 

Then we modify this operator in order to get the further property that for any integrable 
/ E CO((T"\B(xo,2ro)) x M"), one has 



/ Tr{f)dvdx= / fQ{x,v)dvdx. 



(3.36) 



This condition can easily be obtained by considering a regular, compactly supported, nonnegative 
function u with integral 1 in i?(xo,?'o) ^ 1^") cind adding to '7r(/) the function 



fo- f 

T"xR" J(T"\a;)xR" 



U. 
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We obtain a continuous affine operator vr satisfying that for some constant 0,^, one has for any 
integrable / G C^(T^\i3(xo, 2ro)), one has 



||vr(/)|bi <c.||/|bi + 
k(/)||L- <c^||/||l-.+ 



/ - / fodvdx 

(T"\a;)xE" JT"xIR" 



f — fodv dx 



Due to the compact support of u, vr continuously sends L°°((T" \ w) x M"; (1 + |t>|)i'(ix) into 
L°°(T" X M"; (1 + \v\)^ dx), with estimates as above. 

It is convenient to introduce another truncation in time function T such that: 



T = in 



' 100 



and T = 1 in 



48 



(3.37) 



Finahy, we introduce the operator U : C°(([0, T] x [T'^\B{xo, 2ro)] x M^) y qq, r/48] x x 
M^)) ^ C0([0, T] X T2 X M2) given by: 

{Uf ){t, X, v) := (1 - t(t))/(t, X, v) + t(t)[7r/(t, •, OK^;, v). (3.38) 

We finally define V[g] by: 

l^b] := / + n(/|([o,T]x[T2\B(xo,2ro)]xR2)u{[o,r/48]xT2xR2)) in [0)^] X X M^. (3.39) 



3.3 Existence of a fixed point 

The goal of this paragraph is to prove the existence of a fixed point for small values of £, which 
corresponds to the following lemma. 

Lemma 3.1. There exists eo > such that for any < £ < £q, there exists a fixed point ofV 
in Se- 

The proof is almost the same as in [3 Section 3.3]. In order to avoid to repeat it, we only 
give the main arguments and refer to it for the details. We only focus on the main differences. 

We endow the domain with the norm of C°([0,T] X X M2). The existence of a fixed 
point of V on Ss relies on Schauder's theorem. Accordingly, we have to prove that Ss is a convex 
compact subset of C'^([0,T] x x M?), that V is continuous on for this topology, and finally 
that V{Se) C Se- 

That Ss is convex is clear; that it is compact follows from Ascoli's theorem, using both 
uniform Holder estimates and the uniform weighted estimates. 

Now let us discuss the continuity of V. Here the proof of [3 Section 3.3] actually holds 
without further modification. Let us briefly explain the argument. Due to the compactness 
of Ss, it is sufficient to prove that if /„ — )• / in Ss, then V[fn] — V[/] pointwise. Let us fix 
(x, ti) G X M^. Call (X", V^) and {X, V) the characteristics associated to the force F + Vipf'^ 
and F + 'Vip^, respectively. By Gronwall's lemma, {X'^,V"') converges to {X,V) uniformly on 
compacts. 

If there was no absorption (that is, if we took U = 0), then the convergence 

V[fn]{t,X,v)^V[f]{t,X,v) 

would follow from Vip-^"" — )■ Vip^ uniformly on [0, T] x and Gronwall's lemma. The difficulty 
comes from the fact that we have to take into account in V[f]{t,x,v) the various times of 
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absorption on 7". But from the convergence of (X'^,V'^) to {X,V) (uniformly on compacts), 
one can deduce that for n large enough, (X"", 0, x, v) meets 7" the same number of times as 

(X, V){-, 0, X, v), and that the intersection points of {X"", 0, x, v) and 7" converge towards 

those of {X, V){-,0,x,v). Then the continuity of V follows. 

The main point in the proof is to establish that V(5e) C Ss- The crucial estimate here is the 
following. 

Lemma 3.2. Let g G S^, and {X,V) the characteristics associated to F + 'V(p^. Then one has 

\\v\ - \V{t,0,x,v)\ \ < l + t\\F + Vip3\\^. (3.40) 

This lemma is trivial in the case under view, even with \v — V{t, 0, x, v)\ on the left hand side. 
But since the estimate with \v — V{t,0,x,v)\ on the left hand side is not valid in the presence 
of a magnetic field, we prefer to use (j3.40p . 

Let g £ Se- That the point d. is true for V[g] comes from the construction, in particular 
from the choice of the operator 11 (see ()3.36p ). 

Let us explain why the point b. is satisfied by / := V[g]. From the construction, on 7^ one 
has \ f{t'^,x,v)\ < \ f{t~,x,v)\. It follows that 

\fit,x,v)\ < \fo[iX,V){0,t,x,v)]\. 

Now, 

\fit,x,v)\ < \\{l + \v\)yo\\L^[l + \\v\-{\v\-\V{0,t,x,v)\)'^^~'' 

7 

/ — ll/MIt 'riilll X 

< ii(i + bir/oiiL- 

where we used 



1 + \\v\ - \V{0,t,x,v)\\ 
1 + \v\ 



{l + \x-x'\)-^ < ^ + 



1 + 

Note that \\F + V(/7^||oo < ll-^lloo + £ < ll-^lloo + 1- With Lemma l3.2| we deduce that for some 
C > independent of /o and e: 

\{l + \v\rf{t,x,v)\<C\\il + \v\rfo\\Lo^. 
Then the fact that V[g] also satisfies b. follows from the construction of the operator 11. 

Let us now explain the point c. We have the following lemma: 

Lemma 3.3. For g G S^, one has V[g] G C^{Qt\T,t), with := [0, T] x 7". Moreover, for 
any {t,x,v) and (t',x',v') in [0,T] x [T^\a;] x M^, with \v — v'\ < 1, one has, 

\V[g]{t,x,v) - V[g]{t',x',v')\ < C[\\fo\\ci^j2^^2) + ||(1 + W M\ l^. ^^2)] 

X {l + \v\)\{t,x,v)-{t',x',v% (3.41) 

and also 

\V[g]{t,x,v)-V[g]{t,x',v')\<C[\\fo\\ci^j2^^2)^^^^ (3.42) 
the constant C being independent of Jq. 
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This lemma is rather technical. Actually without absorption, this estimate follows from 
Gronwall's lemma and the regularity of V[g] follows from the fact that /q and the characteristics 
are of class C^. But here at each passage in 7", there is a jump between VV[g](t~^ , x,v) and 
\/V[g]{t~ ,x,v). One can see by using an explicit computation based on the last equation in 
(I33I11 that 

I VV[g] {t+,x,v)\<\VV[g]{t-,x,v)\+C\V[g]{t-,x,v)\, 

where V is either V^,. or V^,. 

The main point is that the number n{x, v) of times a characteristic {X, V){t, 0, x, v) can cross 
7~ is estimated as follows. Using dist(7~,7+) > and Lemma [3.21 we infer that 

n{x,v) < C{l + max\V{t,0,x,v)\) < C(l + |t;|). 

This allows to bound VVfs'] using to the uniform estimates on (1 + |v|)'^V[5]. 

Finally, point a. is a consequence of points b., c. and an easy interpolation argument 
between weighted Holder spaces, provided that /o is small enough. 



3.4 A fixed point is relevant 

Let us prove that, provided that e is small enough, the fixed point that we constructed is indeed 
a solution / starting at /o and reaching in \a; at time T. For this we show that y[(7](r) = 
in X R2. 

Call again [X, V) the characteristics associated to F + Vip-^ — S/Tp + £. 
Due to the construction, it is enough to prove the following lemma. 

Lemma 3.4. There exists ei > such that for any < e < ei, all the characteristics {X,V) 
meet 7^" for some time in ^-]. 

Proof of Lemma \3.4\ We denote by {X , V) the characteristics associated to F + £. 

1. We first prove that for all {x,v) G x M?, there exists a G ^] U [^, ^] such that 

X{a,0,x,v) G 7^- := S^{x,v) G 5(xo,ro) x / \v\ > ^ and v.u{x) < • (3-43) 

Let {x,v) G X ^e claim that that there exists t G [f , ^] U ^] such that 

X{t,0,x,v) eB{xo,ro/^), (3.44) 

and 

\V{t,0,x,v))\>a. (3.45) 
We discuss this according to the modulus of V{T/3,0, x,v): 

• If \V(T/3,0,x,v)\ > M > Ml, then one can observe that \v\ > ma,x(m^, 2a), using the 
characteristics equation. Then by Proposition 13. 1| the claim is proved for some t G l"^, 

• If \V{T/3,0,x,v)\ < M, then by Proposition K2[ \V{2T/3,0,x,v)\ > M + 1 > M2, and 
one can once again apply Proposition 13. H to prove the claim for some t G [-^, 
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Now, one can easily see that for some s > with s < ^ < 

X{t,0,x,v) - sV{t,0,x,v) e 5(xo,ro) with V{t,0,x,v).iy < -^\V{t,0,x,v)\, (3.46) 

because a straight Une arising from B{xQ,rQ/2) cuts S{xQ,ro) with angle to the normal ly at the 
circle of value at most 7r/6. The same argument shows that: 

X{t, 0, X, v) - 2sV{t, 0, X, v) ^ B{xo, 3ro/2). 

Now it is clear that, 

\ViT,0,x,v) -V{t,0,x,v)\ < 2s[\\F\\^ + llV^illoo + ||V^3||oo] for t e [t - 2s, t], (3.47) 

|X(r, 0, X, v) - X{t, 0, X, v) + it- T)V{t, 0, x, v)\ 

< 2s^[\\F\\^ + ||V(^i||oo + ||V(^3||oo] for T G [t - 2s, t]. (3.48) 

In the other hand, if C^o is large enough, we have the estimate: 

\X{t - 2s, 0, X, v) - X{t, 0, X, v) + 2sV{t, 0,x,v)\ < y . 

Therefore by the intermediate value theorem that there exists a & [t — ■^,t], such that 
X{t,0,x,v) G S{xQ,rQ). Using (|3.46p . (|3.47|) and (|3.48|) . and provided that Cr^ is large enough 
(in terms of rg only), we deduce that for this a, ()3.43p applies. 

2. Now to prove that all the characteristics meet 7^" during let us compare (X,V) 

and {X, V). Using point a. in the definition of S^, we deduce by Gronwall's lemma and elliptic 
estimates that 

\{X,V)-(X,V)\<Ce. 

Proceeding as previously, we deduce that if e is small enoug h, then for ah (x, v) £ x there 
exists t G ^"Ij such that 

ix,v){t,o,x,v)ej'~. 

□ 

We can now gather all the ingredients to prove Theorem 11.11 

Proof of Theorem Using Lemma 13.11 we deduce the existence of some fixed point / = V [/] . 
Using Lemma |3.4| and (|3.3ip . (|3.32p and (j3.34p . we see that, provided that e is small enough, 
V[/](r) = 0. Hence / satisfies Supp [f{T, ■, •)] C w x M2. 

It remains to prove that / satisfies (jl.ip . This comes from the fact that, due to ()3.12p and 
(|3.34p . one has 

dtf + + (F + Vipf).V,f = in [0, T] x [T'^ \ u] x R". 

Since / is C^, one has 

dtf + v.VJ + (F + V(^^).V,/ = G in [0, T] x T" x 
for some continuous function G. This concludes the proof of Theorem 11.11 □ 
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4 Global controllability for the bounded external field case 

In this section, we prove Theorem 11.21 

We call H a hyperplane in R" such that its image T-L by the canonical surjection s : M"" — t- T" 
is included in lo. We recall that Ti is supposed to be closed. We call nn a unit vector, orthogonal 
to Ti. For / > 0, we denote 

'Hi:='H + [-l,l]7iH- 
Since Ti is closed in T", we can define d G M"^* such that 

C UJ, 

and such that Ad is less than the distance between two different hyperplanes in s~^{'H). 



4.1 Design of the reference solution 

The reference solution is not quite the same as in Section |3l In order to get a global result, as 
explained in Section [21 we will need the following property, refered to as a "non concentration 
property" for the characteristics {X, V) associated to ^ (up to a slight modification inside the 
control zone): there exist c > such that 

yx,y£ T", \X{t,0,x,0) -X{t,0,y,0)\ > c\x - y\. 

The assumption on the control zone uj is motivated by the fact that in this case we can 
atually construct a reference solution whose characteristics satisfy this condition. 



To construct ((/?, /), we start with the following lemma. 
Lemma 4.1. There exists ip G C~(T";M) such that 

A^ = onT'X-Hd, (4.1) 

and 

V<f = nH onT''\nd. (4.2) 

Proof of Lemma \4-1\ In the domain T'^\T-Ld, x i— )• uh coincides with the gradient of a harmonic 
function. Call (p a function in C°°(T"; M), whose gradient coincides in T-L^ with uh] this function 
is automatically harmonic in Tid- Q 



Now given such a if, we can construct if and /. Consider a function y G C^(0, T) satisfying 

SuppJ^C (f,f ), 

y ^ 0' (4.3) 



Set 



/ ^= 


1. 




J[Q,T] 






for t G 


T 






U 


Hy{t)(p{-) for t G 



2T 

— ,T 
3 ' 

T 2T' 
3"' "3" 
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T 



2T 



3 ' 
T 2T 



for t G 0, - U 
Hy{t)nH for t G 

where /x is a positive parameter depending on T and F only, according to the following lemma. 



Lemma 4.2. Given u as above, T > and F, there exists /x > such that all the characteristics 
associated to £ meet 

^3- := {{x,v) G drid X M" / \v\>2 and v.u < -2} , (4.4) 
for some time in where v = inn is the outward unit vector on dTid- 

Once defined 7p, we define / : [0, T] x x as previously by ([322]) -([323]). 
Proof of Lemma \4-S\ Let (x, G T" x M". Call (X, V) the characteristics associated to £. We 

•T 
. 6 ' 

/T 



discuss according to 1^(^,0, • nn- 

• If y(^,0, • njy is large enough, say larger than c > 0, then one sees easily using the 



characteristic equation that there exists t G [g-, such that {X ,V){t,Q,x,v) G 7 . 

For the other {x, v), one can find > such that 0, x, v) ■ nn > c. Then there exists 

i G [^,^] such that (X,V){t,0,x,v) G 7^". 

□ 



4.2 Definition of the fixed-point operator 



For A G (0, 1], we define again a subset of C^\Qt) on which we will define the operator V 
(which actually depends on A): 

a- II /Mn(5 -7)c?i'llc*i(nT) - ^' 

b- 11(1 + |u|)'^(g -7)||loo(q^) < Ci ||/o^||ci(T"xR") + 11(1 + kl)'^/o'llcO(T"xR") 

C- lb - l\\ch.(^Q^-^ < C2 ||/o^|lci(T"xM") + 11(1 + \v\T'fo\\cO{T"xR") 

d. VtG[0,T], /^„^^„5(t,x,i;)dxdi; = /[p„^^„/o^(x,'y)(ix(iu I , 



(4.5) 



with ci, C2 to be fixed later depending only on 7, T and a; (and hence on (/, if)), but not on A; 
here, di and 82 are fixed as follows 



Si 



7 — n 



2(n + l)(7 + l) 



and 82 



7 + 1 



For fixed ci and C2 large enough depending only on (/, ip), one has for A small enough depending 
on e that 



f^dvdx 



see (j2.3p . Hence in that case g{t,x,v) = /q(x,u) + f{t,x,v) belongs to for A < /u(e), so 
7^ 0. From now, we suppose that this is the case. 
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We write Ti := Ti — dnn, T2 := T-L + dnu and F := Fi U F2. Let v = —Uh on Fi and u = rih 
on F2. We define 



7 

,2- 



:= G F X / < -1} , (4.6) 

:= |(2;,7;) G F X M" / \v\ > 1 and vm{x) < -3/2| , (4.7) 

7+ := |(x,u) G F X M" / 7;.i/(x) > 0} . (4.8) 

Note that 7^" defined in ()4.4p can be reformulated as 

^{x,v) G F X ]R" / |t;| > 2 and vm{x) < -2| . 



Again, we observe that 

dist((F X M")\7-;72-) > 0. 

We introduce a C°° H regular function [/ from F x R" to M the same way as previously, by 

0<U <1, 

f7 = 1 in (F X M")\7-, (4.9) 
U = 0m 7^-. 

The function T is again introduced by (|3.32p . As in Section [3l we define vr as a continuous affine 
extension operator vf from 0^(712(1]^) to C''(T"';R), and which has the same property that 
each C"^-regular function is continuously mapped to a C"-regular function, for any a G [0,1]. 
Moreover, we manage again in order that for any / G C^{7i2d',^)j (|3.36p occurs. The operator 
n is given by ()3.38p . 

Now, given 51 G 5^, we first define ip^ by ()3.33p . 

Then we introduce f = V[g] as the solution of the following system: 

' f{0,x,v) = on T" X M", 

dtf + v.V.f + {F^ + V{ip9 - if) + ^iy{t)nH).V,f = in [0,T] x [(¥" x M")\7-], (4.10) 
[ f{t,x,v) = [l-T{t)]f{t-,x,v) + T{t)U{x,v)f{t-,x,v) on [0,r] X7-. 

The meaning of this equation is the same one as in Section [3] (and fJ.y{t)nH plays the same role 
as £ in Section [3]). Recall that was defined in (j2.5p . 
Then, as for Section [3l we define V[g] by 

V[g] :=7 + n(/|[o,T]xW2,xM"u[o,T/48]xT"xR") m [0,r] X T'^ X M". (4.11) 

Again, /|[o,T]x'H2dxK"u[o,T/48]xT"xM" is regular, and, together with the construction of FT, it 
win follow that V[g] is in C^{[0,T] x T" x M"). 

Considering the form of (|4.10p . the characteristics that we consider in the sequel are (X^ , V^) 
associated to F^ + V{(p^ — ^) + fiy{t)nH, which coincide with the ones associated to F'^ + V(p^ 
outside the control zone, but not necessarily inside. 
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4.3 Existence of a fixed point 

Now our goal is to prove the following lemma. 

Lemma 4.3. For any small e > 0, there exists \{e) > such that for any positive A < A(e), 
the operator V has a fixed point in . 

Proof of Lemma \4-'J\ We prove Lemma 14.31 bv checking the assumptions for Schauder's fixed 
point Theorem on V. We will sometimes forget the indices and exponents e and A. 

1. Again, 5 is a convex compact subset of C^{Qt)- 

2. The continuity of V can be proven in the same way as in Section 3. 

3. The difficulty is to check that for A small, one has V{S^) C S^. Accordingly, we have to 
check the points a., b., c. and d. for V[g\. 

That V[g\ satisfies d. comes directly from the construction. That V[g\ and consequently V[g\ 
satisfies estimates as b. is not difficult and proven as in Section [3l In particular Lemma 13.21 is 
still satisfied. 

For what concerns point c. we have as previously (see also [3 Lemma 4, p. 370]) 

Lemma 4.4. For g G , one has V[g\ G C^{Qt\^t), with := [0,T] x 7". Moreover, for 
any {t,x,v) and {t',x',v') in [0,T] x [T^\a;] x M^, with \v — v'\ < 1, one has, 

\V[g]{t,x,v) - V[g]{t',x',v')\ < C[||/o||c1(t2xr2) + ||(1 + bir+VollL^CT^xR^)] 

X {l + \v\)\{t,x,v) - it',x',v% (4.12) 

and also 

\V[g]{t,x,v) - V[g]{t,x',v')\ < C[||/o||c1(t2xR2) + 11(1 + H)^^^ fo\\L--{T2,^K2)]\{x,v) - {x',v')\, 

(4.13) 

the constant C being independent from /q. 

The central part is point a., where the smallness of A and the non concentration property of 
Tp are used. We begin by a lemma which asserts that the non concentration property is preserved 
by a small perturbation. Recall that {X^ , V^) are associated to F"^ + V((^^ —'ip)+ ^y{t)nH- 

Lemma 4.5. There exists c > such that for any A small enough (in terms of T , co and F), 
for any g £ S^, one has 

V(x, y) G (T")2, yt G [0, T], c" V - y\ < \X3{t, 0, x, 0) - X^it, 0, y,0)\ < c\x - y\. (4.14) 

Proof of Lemma \4-5\ Define {X, V) as the characteristics associated to the force iiy{t)nH- It is 
clear that {X,V) satisfy the non concentration property: 

V(x,y) G (T")2, Vt G [0,r], \X{t,0,x,0)-X{t,0,y,0)\>\x-y\. (4.15) 

(This is actually an equality!) Now, it follows from Gronwall's inequality that for a constant C 
depending only on 3^ and F, one has 

\\{X^,V^) - (X,F)|bo([o,T]2xT"xR'^) < C(e + A2). (4.16) 
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One can get a further inequality in the following way (when it is not explicit, the norm considered 
is the L°° one) 

^ II V(X3, V9){t, s, X, v) - V(X, F)(t, s, X, v)\\ 
< \\VV3{t,s,x,v)-VV{t,s,x,v)\\ 

+ \\V.^Eg{t, X3{t, s, X, v))VX3{t, s, X, v) - V^.Sj(t, X{t, s, x, v))VX{t, s, x, v)\\ 
+ II V,,,F^(t, X^it, s, X, v), Va{t, s, X, v))V{X3,V3){t, s, x, v) 
- Vx,vF^{t, X{t, s, X, v),V{t, s, X, v))V(X, V){t, s, x, v)\\ 

where V stands either for V^; or for V^,. Now the second term is bounded as follows 

\\VxEg{t,X3{t,s,x,v))VX3{t,s,x,v) -VxEj{t,l{{t,s,x,v))VX{t,s,x,v)\\ < Ai + A2 + A3, 

with 

Ai = \\V^Eg{t, X9{t, s, X, v))VX9{t, s, X, v) - V^Eg{t, X9{t, s, x, v))VX{t, s, x, v)\\, 
A2 = \\VxEg{t,X3{t,s,x,v))\7X{t,s,x,v) - \7r,Ej{t,X3{t,s,x,v))\7X{t,s,x,v)\\, 
A3 = \\VxEj{t, X9{t, s, X, v))VXit, s, X, v) - V^Ejit, X{t, s, x, v))VX{t, s, x, v)\\. 

Now 

Ai < \\VccEg\\co(^n^)\\VX3{t,s,x,v) - VX(t, s, x, i;)||c.o([o,r]2xT"xR"), 

A2 < llVx^'g - Vx-i^j||c'0(nj,)||VX||(^0([o,r]2xT"xR")5 

A3 = 0. 
Hence we obtain 

\\VxEg{t, X9{t, s, X, v))VX^t, s, X, v) - VxEj{t,X{t, s, x, v))VX{t, s, x, v)\\ 

< C{e+\\VX3{t,s,x,v) - VX(^,s,x,^;)||c.O([o,T]2xT"xK"))• 
We treat the term concerning in the same way and obtain 

II V.,.F^(t, X3{t, s, X, v), V9{t, s, X, v))V{X9, V9){t, s, x, v) 

- V^,^F^(t, X{t, s, X, v),V{t, s, X, v))V(X, V){t, s, x, v)\\ 

< C{X + ||V(X^ V9){t, S, X, V) - V(X,V)it, S, X, l^)lbo([0,T]2xT'^xKn))- 

It follows then by Gronwall's lemma that for a certain constant C, one has 

\\{X9,V9) - (^,"V^)||loo([o,t].ci(T"xM")) < C{e + A). 

Hence, if e and A are small enough, then (j4.15p is still valid when replacing (X, V) by {X^ , V^), 
up to a multiplicative constant. This gives ()4.14p . □ 

Let us come back to the proof of point a. Let us treat the L°°-norm; the C^^ one will follow 
by interpolation. From (jHH), we deduce that X9{t,0,-,0) : T" T" is invertible; call {X[)-^ 
its inverse, and define the function Wf : [0, T] x T" — )■ M" by 

Wfix) := yf(t,0,(Xf)-Hx),0). 
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One can describe {X^)~^{x) as the initial position of a particle, which starting with velocity 
0, reaches x at time t; then Wf{x) is its velocity at time t. 
Let us give an estimate on v — Wf{x). First, 

v-Wf{x) = V3{0,t,X3{t,0,x,v),V3{t,0,x,v)) - V'^{t,0,{X^)-^{x),0). 

By using Gronwall's lemma on V{0,t,-,-), we deduce that for some constant independent of 
A G (0, 1] 

\v - (x)| < C {\X9{t,0,x,v) - {X^)~\x)\ + \V9{t,0,x,v)\) . 

That the constant is independent of A comes from the fact that we have uniform Lipschitz 
estimates on + V((/?^ - ^) + tJCP{t)nH for A G (0, 1]. 

To estimate the first term, we first notice that the non-concentration property ()4.14p gives 

{clr\Xlr\x)-Xa{^,t,x,v)\ < |X^(t,0,(Xf)-i(a;),0)-X^(t,0,X5(0,t,x,i;),0)| 

= |x - X'^(t,0,X^'(0,t,x,7;),0)| 

= \X3{t,Q,X3{Q,t,x,v),V3{{),t,x,v)) - X3{t,{),X3{id,t,x,v), 

where the first equality comes from the definition of {Xf)~^, and the second one of the flow 
property. 

Now Gronwall's lemma for X^(t, 0, •, •), we deduce that for some constant C > independent 
of A G (0, 1] one has 

\Xf^{t,0,X>^{0,t,x,v),V3{0,t,x,v)) - X3{t,0,X^{0,t,x,v),0)\ < C\V'^{0,t,x,v)\. 

Finally we deduce that for some constant K > independent of A, one has, for any A G (0, 1] 
and any g ^ S^, 

\v-Wt^{x)\ < K\V3{0,t,x,v)\. (4.17) 



Now, one has 

\f{t,x,v)\ < \f^[{X3,V<^){0,t,x,v)]\ 

< A2-"||/o(1 + |^;|)^||ioo(Trn,M") (l + ^ \V^{0,t,x,v)\^ 

Using ()4.17p . we get that 

1 



■7 



\f{t,X,v)\ < A2-"||/o(1 + ( 1 + - Wt{x)\ 



It follows that 



/ fit,x,v)dv\<X''-mo{l + \v\r\\L^^jn^^„J' (l + ^\v-Wf{x)\) \v. 



We deduce that 



V[g]{t,x,v)dv\ < kA^-"||/o(1 + |z;|)^||L^(TnxR'^)i^"A". 



One deduces from the construction of V that 

V 2-n 



{V[g] - f){t,x,v)dvh^^nr) < CA^""ll/o(l + \v\r\\L^iTn^u^~,\^ < C(/o)A^. (4.18) 
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Now we turn to the Holder estimate. It follows by interpolation between points b and c, that 
for a certain constant C independent from A, and for 7 = ^^^^-j^ and 5 = 7/(7 + 1) one has 

\V{g\ -1\]<C [||/o^|bi(T"xM'^) + 11(1 + bir/o I|C0(T"XM")" • 
We deduce that, for A < 1 and another constant C (depending on /q but not on A), 

{m-l)dv\\c^(^n^^<C\^-r 

Now we interpolate again this inequality with (|4.18p . We get that for 5i, one has 

II j{y[g\-l)dv\\cs,^^^^<K'\. 

which concludes the point a, for it is sufficient to find a proper A. This finally proves V{S^) C S^. 

□ 

4.4 A fixed point is relevant 

Now we can prove that the characteristics associated to the fixed point are relevant: 

Lemma 4.6. There exists ei > such that for any < e < ei, all the characteristics [X,V) 
meet 7 for some time m [i^g, ^^J- 

Proof of Lemma We recah that by the scaling = \'^F{Xt,x, j), so that ||-^'^||Lf°^ ^ < 
^^ll-^IUj"^ „ ■ As for Lemma 13.41 the proof follows, recalling the Gronwall's estimate (|4.16p . and 
the fact that the characteristics associated to the reference solution / meet 7^" x (We 
recall that /i was defined when we have constructed the reference solution /.) □ 

Finally, we can conclude the proof of the theorem. 

Proof of Theorem Using Lemma 14.31 we deduce the existence of some fixed point g for 
A sufficiently small. Using Lemma 14. 6| and (|3.32p . ()4.9p and ()4.10p . we see that it satisfies 
Supp [g{T, •, •)] C a; X . Now, we define /(t, x, v) = g{j,x, \v), which satisfies the conclusions 
of Theorem 11.21 The fact that (jl.ip is satisfied for some G supported in iv is done as in Section 
El □ 

5 External magnetic field case 

In this section, we prove Theorem 11.31 that is the local controllability result for the external 
magnetic field case. 

5.1 Rephrasing the geometric assumption 

We begin by transforming the geometric assumption (jl.lip in a way that is easier to handle in 
the sequel. For K a compact subset of and r > we denote 

Kr := {x G / d{x, K) < r}. (5.1) 

The geometric assumption can be reinterpreted with the help of the folllowing lemma. 
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Lemma 5.1. Let K CT'^ such that b > on K and satisfying (jl.lip . Then there exists b>0, 
d > and D > such that 

b>b onK2d, (5.2) 



Vx e T^ Ve e S\ 3te [0, D], Vs g 



d 

t,t + - 
' 2 



X + se £ Kd- 



(5.3) 



Proof of Lemma \5.1\ An easy argument relying on the compactness of shows that for d > 
suitably small, one has (|5.2|) . 

To prove ()5.3p . we use the compactness of x S"*^. For any {x,e) G x S^, there exists 
t G M"*" such that x + te £ K. One deduces that for (x', e') in an open neighborhood of (x, e) in 

X §\ one has x' + te' G /sTrf/s. 

Hence by compactness of x S^, there exists a maximal time D such that for any (x, e) G 
T'^ X there exists t G [0, D] for which x + te G Kd/2- Now if x + te G and x + t'e ^ i^^, 

then one has |t — t'| > d/2, since dist(if^/27 \ -^d) ^ '^/2. The conclusion (|5.3p follows. □ 



5.2 Design of the reference solution 

The first step consists in building the reference solution, once again distinguishing between high 
and low velocities. We first treat the case of large velocities. We prove that with the geometric 
assumption on b, high velocity particles spontaneously reach the arbitrary open set. One can 
observe that this is very different to the case of bounded force fields. Actually we can prove a 
stronger result than announced, since we can add to the Lorentz force any additional bounded 
force field. Such a generalization will be actually crucial for the proof of Lemma 15.31 

Proposition 5.1. Let T > and ro > 0. Let b satisfy the geometric condition (jl.lip . 
There exists m G M~^* large enough depending only on b, T and ui such that for all 5 G 
L°°{0,T; tyi'°°(T2xM2)) satisfying \\d\\L°- < 1, the characteristics (X,V) associated to b{x)v-^ + 
5" satisfy: 

Vx G T^Vt; G such that \v\ > m, 3t G (r/4,3r/4), X(t,0,x,t;) G B{xo,ro/2) 

and for all s e [0,T], y < |F(s, 0, x, i;)| < 2|i;|. (5.4) 

Proof of Proposition \5.1[ We prove Proposition 15. II in several cases of increasing complexity. In 
a first time (Cases 1-3), we suppose that 5 = 0. In Case 4, we explain how to take J into 
account. 

In all cases, we define 

b := max6(x). (5.5) 

XGT2 

1. An enlightening case: constant magnetic field modulus. Let us first suppose b constant; for 
readability we assume here that 6(x) := 1. 

As noticed in [T] Appendix A, p. 373-374], there are only a finite number of direction in 
§1 (identifying with [0,27r[, we denote them Q:i,...,a7v £ [Oj^ttO for which there exists a 
half-line in which does not intersect B{xo,rQ/8). Indeed if the slope is irrational, then each 
corresponding half-line is dense in the torus, and consequently meets i?(xo,ro/8). If the slope 
is rational, say p/q with p G Z, g G N \ {0} and gcd(p, q) = 1, then these half-lines L are closed 
periodic lines in T^. Due to Bezout's theorem, the distance between to consecutive lines in 
s~'^{L) is less than min(ji| , |), and the conclusion follows. 
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We introduce the neighborhoods of a^: 

Vi = {a,-Pi/2,ai + Pi/2), 
as fohows. Let /3j > and t <T small enough such that 



T T 

Pi < — and — < min d{Vi,Vj] 
4 4 ij^j 



By a compactness argument, there exists a length L > such that for any x S T^, Va^ E 
Vi, any particle starting from x with a direction has to travel at most a distance L 
to meet B{xQ,ro/8). 

We fix m large enough such that: 

Tm:=-< t/4.. 
m 

This is the time "free" particles with velocity m take to cover the distance L. We observe that 
for any \v\ > m, we have T|^j := -j^ < T^- 

Now let X S T^,u S with \v\ > m. Let us discuss according to the direction of v. 

• First case : |^ G U^^^ Vj. 

We denote (X*, V^) the characteristics associated to free transport. 
We have, for any t < Tj^j, 

\X*{t + T/4,T/A,x,v) - X{t + T/A,T/4,x,v)\ < \v\^ = < 

2 2\v\ 2m 

We can impose m large enough such that ^ < '''o/S- As a result: 

3t G (r/4,r/2],X(t,0,x,t;) E 5(xo,ro/4), 
and (|5.4p is trivial here since 0, x, is conserved. 

• Second case : G U^^^Vj, say Vj. 

The idea is to simply wait for a time r/4. Let us consider 

{x',v') := (X{{T + T)/A,T/4,x,v),Vi{T + T)/A,T/4,x,v)). 

We observe that because of the "rotation" induced by the magnetic field and due to the 
choice of f3i , 

and thus we are in the same case as before. 
Consequently we have proven that: 

3t e iT/4,3T/i],X{t,0,x,v) £ B{xo,ro/i). 

2. Positive magnetic field modulus. Here we suppose that 6 > on T^. 
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We are in the case where in Lemma [5.1[ we can take K = = and 

b = inf h. 

Keeping the same notations as before, we set r G (0, T] and /3j > in order that 

Pi<b^ < mind(Vi,Vj). 
The proof is very similar to the previous one. Indeed, the following estimate still holds: 

\X*{t, r/4, X, v) - X{t, r/4, x,v)\< ^b. (5.6) 

Let X G T'^,v € M^. We distinguish as before between two possibilities. Using the previous 
inequality ()5.6p . the first case holds identically for m large. For the second case just have to 
check that with this magnetic field, the velocity is rotated by an angle at least equal to /3j after 
some time t G (0, j). 

We use the following computation for general {x,v). Denote by 9{t) the angle (modulo 27r) 
between and V{t,0,x,v). Taking the scalar product with V{t,0, x,v) in: 

dV(t,0,x,v) , NN-TT/ ni 
^ ' = b{X{t, 0, X, v))V{t, 0, X, v)^, 

we obtain that |y(t, 0, x, v)| = \v\. Then, taking the scalar product with , we obtain: 

s\ue{t)e'{t) = b(X{t,0,x,v))sm6{t), 

so that 

d'{t) = b(X{t,0,x,v)), (5.7) 

(even if sm9{t) = in which case one considers the scalar product with v.) We deduce that 
d'{t) > b. 

Thus going back to {x,v), by the intermediate value theorem and the definition of the 
neighborhoods Vj, there is a time Tq less or equal to t/4 for which we have: 

F(To + |,|,x,7}) GSi\uiIiH, 

and we conclude as previously. 

3. Magnetic field modulus satisfying the geometric condition. Let us consider the general case 
for b, but without the additional force ^. 

Given K satisfying the geometric condition we introduce d and D as in Lemma l5.ll 

Let 

[/ := T2 \ Kd, 

where we recall the notation (|5.ip . We assume here that r E (0, T] and /3j are such that 

We denote by (X*,!/*) the characteristics associated to free transport, while {X,V) corre- 
sponds to those associated to the magnetic field. 
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Let X G T2,w G We once again distinguish between the two possibihties. As before the 
first case is still similar since ()5.6p is still valid. We have to give a new argument for the second 
case. 

We will assume that m is large enough so that Tm < §• We distinguish between several 
sub-cases: 

a. Assume that X{t,0,x,v) G for some t in a time interval of length at least equal to ^ 
inside [j, Then one can apply the positive magnetic modulus case (case 2). 

b. Assume more generally that C^{{t G [^,^],X{t,0,x,v) G Ka}) > T/4. On U, one has 
6 > 0, so the angle of V{t, 0, x, v) with v is non decreasing over time. It follows that we 
can apply (|5.7p to each passage of the particle in and we conclude as before. 

c. We assume now that the previous cases do not hold. Then X{t,0, x,v) remains in T'^\K(i 
at least during a time ^ in (j, ^). 

By ()5.3p . each passage in T^\Kd of X*{t,0,x,v) lasts at most D/\v\. Actually, in U, 
the characteristics X are not straight lines since they are modified by the magnetic field. 
Let us prove nevertheless that if \v\ is large enough, then the particle can remain at most 
during a time D/\v\ in U. 

Let X G [/, and G Si, let o" G ^). By Lemma l5.ll there exists < -[^ such that 
A'*(cr + s, a, x, v) G K. Now we can evaluate as for a previous computation: 

X*(o" + s, a, x, v) — X{a + s, a, x, v) 

We can choose m large enough such that for any > m, X{a + s,a,x,v) G K^- Hence 
at each passage of X{t,0,x,v) in T'^\Kd lasts at most during a time D/\v\, which proves 
the claim. 

This involves that there are at least [-^^J — 1 passages in U, and therefore there are also 
at least [^^J — 2 passages in K^. This is larger than for \v\ large enough. 

Now we denote by t' a time for which X{t', 0, x, v) G K^i, with X{t, 0, x, v) ^ for t < t' 
and t close to t' . Let us show that X(t' + s,0,x,v) remains in for s < if the 

velocity is large enough. We have for all s G [0, ^y^], 

\X*{t' + s,t',x,v) -X{t' + s,t',x,v)\ <b\v\ ^^^^ . 

On the other hand, by (|5.3p . each passage of X* in lasts at least Hence we can 
choose m large enough such that for any \v\ > m, X{t' + s,t',x,v) G for s G [0, 

Consequently, X{t, 0, x, v) remains in during a time inside ( j, ^), and we conclude 
as before. 




4. With a nontrivial additional force 5- 

Let us finally explain how one can take ^ into account. First, we consider the equations for 
\V\ and 9, where 6 is the angle between v and V{t, 0, x, v). The following computations are valid 
for V large so that |y(f,0, does not vanish and for a time interval where 6 G [— 7r/2, 7r/2]. 
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• For what concerns \ V\, it suffices to take the scalar product with V{t, 0, x, v) of the equation 
of V . We infer 

^|F(t,0,x,^;)|2 =25-F(t,0,x,t;), 
at 



so that 



In particular, for m large enough, one has for all {x,v) G x with \v\ > m, 

^ < \V{t,0,x,v)\ <2\v\. 



(5.8) 



(5.9) 



For what concerns 6, taking the scalar product of the equation of V with v we deduce 
d 



Hence 



\V{t,0,x,v)\ \v\cose{t) -\V{t,0,x,v)\\v\e'{t)sme{t) 



b{X(t,0,x,v))V {t,0,x,v) -v + d -v. 



\V{t,0,x,v)\\v\9'{t)sme{t) 



= b{X{t, 0, X, v))\V{t, 0, X, v)\\v\ smi9{t)) -d-(v- ^M'^'^)H cos Bit) 

V \V{t,0,x,v)\ 

We notice that 

V{t,0,x,v)\v\ V{t,0,x,v) „ . . . 

^ - imo,.,.)i = ' - fm:^ ^""''^ = pm.o..)}-(-)' 

where P|y(f q x i;)}^ (^) denotes the orthogonal projection of v on {^(t, 0, x, f )}^. So 



?'(t) = 6(X(t,0,x,t;)) + 



1 



P{V{t,0,x,v)}^ 

\V{t,0,x,v)\" |i;|sin6'(t) 



7^ 



(5.10) 



Note that 
so that: 



\P{Vit,o,x,v)}^i'")\ = 1^1 |sin(0(t))|, 



1 



\V{t,0,x,v)\ 



P{V{t,0,x,v)}^('") 



\v\ sin 6{t) 



1 



\Vit,0,x,v)\ 
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Now let us revisit the three sub-cases of Case 3 to include ^. 



a. Assume that X{t,0, x,v) G for all t in a time interval of length at least equal to 
Then using (|5.9p and (|5.10p we deduce 



9'{t) > 6-2 



m 



(5.11) 



so one can conclude as in the positive magnetic modulus case. 
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b. Assume more generally that C^{{t £ [^,^],X{t,0,x,v) G Ka}) > T/4. On U, one has 
5 > 0, so the angle of V{t, 0, x, v) with v satisfies 

e'it) > --lisiloo, (5.12) 

m 

and (|5.1ip when X{t, 0, x, v) € K^- In total the variation of 6 is no less than ^ — ^||5^||oo) 
so one can conclude as previously (taking m large enough). 

c. We assume now that the previous cases do not hold. Then X{t, 0, x, v) remains in T^\i^d at 
least during a time ^ inside ( ^). Let us compare the characteristics {X, V) associated 
to ^+b(x)v'^ with the characteristics {X, V) associated to the magnetic field b{x)v^ alone. 

Let x €z U, and G Si, and let let cr G (-j, ^). Using the analysis of case 3, there exists 

t' < such that X{a + t',a,x,v) G K^. Now comparing {X,V) and {X,V) and using 
Gronwall's inequality we deduce 

\V{a + t',a,x,v) -V{a + t',a,x,v)\ < ||5||oo exp(||6||;yi,co (1 + 2|t;|)t'), 
\X{a + t',a,x,v) - X{a + t',a,x,v)\ <t'\\d\\oceM\Mw^-^i^ + '^M)t')- ' 

Using that \v\t' is of order 1 and taking m large enough, we see that for any |f | > m, 
X{a + t',a,x,v) G i^3d/2. Hence each passage of X{t,0, x,v) in T^\X3^/2 lasts at most 
D/\v\. We deduce as previously that there are at least [^^J — 2 passages of X{t,0, x,v) 
in K3a/2 during (j, ^). 

Now reasoning as in Case 3, using Gronwall's estimate (|5.13p . we see that if X{a, 0, x, v) G 
I^3d/2i and m is large enough, then X{a + 1', 0, x, v) remains in K2d for all times t' < 
and we conclude as before. 

□ 



Now let us turn to the case of low velocities. This time we proceed as in the case of bounded 
force fields and prove that an analogue of Proposition 13.21 holds: 

Proposition 5.2. Let t > and M > 0. There exists M > 0, £ £ C~([0,r] x T'^;R'^) and 
(f G C°°([0,r] X T2;M) satisfying 

£ = -Vip m [0, r] x (T^Bixo, ro)), (5.14) 
Supp{£) C {0,t) xT"^, (5.15) 
Aip = m [0, r] x {T\B{xo,ro)), (5.16) 

such that, for any ? G L~(0, T; VF1'°°(T2 x R^)) satisfying < 1, «/ {X,V) are the char- 

acteristics corresponding the force ^ + £ -\- b{x)v-^, 

y{x, u) G X B{0, M), V{t, 0, x, v) G B{0, M) \ B{0, M + 1). (5.17) 

Proof of Proposition \5.4A Again, we introduce 9 and £ as in the proof of Proposition 13.21 Again, 
one can choose C and then r' such that 

V(x, v)eT^ X B{0, M), V{t, 0, x, v)eR^\ B{0, M + 2 + r||5||oo). 

Let us denote by {X, V) the characteristics corresponding to the force £ alone. 
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We first observe that we have: 

j^\V\'' = {^{s,X,V)+£{s,X)).V. 
Thus, using Cauchy-Schwarz and Gronwall's estimates, we obtain: 

|y|2<max(l,|^|2e*(M-+ll^ll-)). 

We evaluate: 

\X{t,{),x,v) -X{t,Q,x,v)\ < [ \V{s,0,x,v) -V{s,0,x,v)\ds, 

Jo 

\V{t,0,x,v) -V{t,0,x,v)\ < [ \\£{s,X{s,0,x,v)) - £{s,X{s,0,x,v))\ 







+ \mX,V)\+b\V{s,0,x,v) 



ds 







< / \\^£\\oo{t- s)\V{s,Q,x,v) -V{s,Q,x,v)\ds 

2M 

l^^lloo + ll^ll 



+ max ( T, ^^11 (el(ll^ll°°+lli^ll°°) 



By Gronwall's inequality: 



\V{t, 0, X, v) - V{t, 0, x,v)\< max [T/2, ^^'^ (e^(l|g|l°°+llg|l°°) _ i)\ e^W^'^W 

V ll'-'lloo ~l~ llulloo / 

(5.18) 

For t = t' , we have: 

r' r' 

where C and C depend only on uj,M, and the conclusion follows as previously since 

I \V{t,0,x,v)\ - \V{t',0,x,v)\ I < |r-r'|||5|U. 

□ 

The reference solution. Let us now describe the reference solution. Consider xq in w and 
ro > such that B{xq, 2ro) C OJ. We define the reference potential Ip : [0, T] x — t- M as follows. 
We apply Proposition 15. II with r = T/S, we obtain some m > such that (|5.4p is satisfied. Then 
we apply Proposition 15.21 with r = T/S and 

M = max (m + 1^,100, 32ro(& + 1)) , (5.19) 

and obtain some lp2, £2 and some M > such that (j5.17p is satisfied. We set 

w.-) -- 



and 



£{t,.) 



forte [0, f]u [f^,r], 

^2(t-i,-) fort G [f,f ] 

forte [0, f]u [f^,r], 



£:2(t-^,-)fort€[^,f]. 
Then once defined ^, we define / : [0,r] x x as previously by 
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5.3 Proof of Theorem [TTBl 

We consider Ss the same convex set as in the proof of Theorem 11.11 and V the same fixed point 
operator with F = b{x)v^. As before, the proof consists in proving first the existence of a fixed 
point, and in a second time in proving that such a fixed point is relevant. 
For what concerns the existence of a fixed point we have: 

Lemma 5.2. There exists Eq > such that for any < e < Eq, there exists a fixed point ofV 
in Ss- 

Proof of Lemma \5. M The proof of Lemma [5. 21 is exactly the same as the one of Lemma [3. II and 
is therefore omited. Note in particular that a variant of the crucial Lemma [3.21 is still valid here, 
using dSSl). □ 



In the second part of the proof we show that a fixed point is relevant. In this part lies the 
main difference with Theorem 11.11 This is given by the following lemma. 

Lemma 5.3. There exists ei > such that for any < E < Ei, all the characteristics {X,V) 
associated to b{x)v-^ + £ — + V^p^, where f is a fixed point of V in S^, meet 7^" for some 

r rr 1 1 T'-i 

time m [j^, -jyI- 
Proof of Lemma 15. At 

We begin by noticing that \/(p^ — VTp satisfies 

WVif^ - V^lloo < 1, (5.20) 

provided that e is small enough, which we suppose from now. Consequently we can apply 
Propositions 15. ll and 15.21 to 5 := '^'^^ — S/Tp. 

It follows that any (x, i)) G x is (at least) in one of the following situations: 

• If |y(^, 0, x,v)\ > M, then using (|5.20p . we deduce \v\ > m. Hence there exists r G ^] 
such that 

X(r,0,x,7;) GB(xo,ro/2), (5.21) 
and reasoning as for (j5.9p we deduce that for all s G [0, -1^] one has 

M 

\V{s,Q,x,v))\>—, (5.22) 

where M was defined in (|5.19p . 

• Or \V(^,^,x,v)\ < M, so \V{^,0,x,v)\ > M + 1, and there exists r G [j^,^^] such 
that (f572T]l is true and ([OS]) is valid for ah s £ [^,T]. 

Let us consider (x, v) in the first situation, the reasoning being identical for the second situation. 
As in the proof of Lemma I3.4[ we deduce the existence of some s > with s < ^ < ^ , 

X{t,0,x,v) - sV{t,0,x,v) G S{xo, ^) with V{t,0,x,v).i' < -^\V{t,0,x,v)\. (5.23) 

Let us show that this involves for \v\ large enough the existence of r* G [r, t] such that 
X* := X{t, 0, X, v) — (r* - t)V{t, 0, x, v) G S{xo, tq). 
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We have for a G [r — s, r]: 



M 

— < \V{(T,Q,x,v)\ < 2\v\ 



V{a,Q,x,v) y(r,0,x,w) 



\V{a,Q,x,v)\ |y(r,0,x,^;)| 



< s 



h + 



M 



\X{a, 0, X, v) - X(r, 0, x, v) + {t - a)V{T, 0, x,v)\< —i2\v\ + || 



/I 



(5.24) 
(5.25) 
(5.26) 



Estimate ()5.25p comes from the identity 



d ( V{a^ 0, X, v) 



da \\V{o,^,x,v)\) \V{p,^,x,v)\ 



^{a,^,x,v) V^^(a,x,t;)-y(a,0,x,7;) 



+ 



|F(cT,0,x,t;)|3 



1/(0", 0, X, u) 



Let us check that this involves the existence of t S [r, r — s] such that (X(t, 0, x, f ), F (t, 0, x, u)) G 
7^~. The existence of of t G [r, r — s\ such that X(t,0, x,t>) G S'(xo,?'o) fohows from (|5.26p and 

•s^, , , „^ f„ s 8ro , , , , 2M + 1 24ro rn 

-(2|^;| + ||V(^^||oo) < ^(2^1 + 1) < 8ro _ _T < — ° < ^. 



M2 



M 



At such a t, from (fO¥|) . we have \V{t, 0, x, ?;)| > 2 since M > 4. 

The fact that at such a moment t, one has y(t, 0, x, u ).i/(X(t, 0, x, f )) < — ■g|y(t, 0, x, )| 
comes from 



F(t,0,x,t>) ^(r,0,x,t;) 

• v\X(t, 0, X, u)) - — — — -j- • z^(x*) 



|F(t,0,x,^;)| 



< 



|y(r,0,x,t;)| 
y(t,0,x,i;) 



y(T,0,x,7;) 



|F(t,0,x,T;)| |y(r,0,x,7;)| 



+ 



i/(X(t, 0, X, f)) — i^(x* 



^ (^+lF)tT + -|^(*'0'^'^)-^* 



< 



,4ro 24 1 



and from (|5.23p . This concludes the proof of Lemma 15.3 



□ 



Let us finally gather all the pieces to prove Theorem 11.3 



Proof of Theorem \l.'J[ Using Lemma [5.2[ we deduce the existence of some fixed point / of V in 
Se- Using Lemma 15.31 we can again use the definitions ()3.3ip . (|3.32p and (|3.34p to deduce that 
Supp [fir, •, •)] C X and one checks that / satisfies the equation for some G as previously. 
This concludes the proof of Theorem 11.31 

□ 
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